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Abstract 

rn"; 

It is known that the Swift-Hohenberg equation du/dt = -(d 2 + l) 2 u+s(u-u 3 ) can be re- 
duced to the Ginzburg-Landau equation (amplitude equation) dA/dt = 4d 2 x A+e(A-3A\A\ 2 ) 
by means of the singular perturbation method. This means that if s > is sufficiently 
Qh. small, a solution of the latter equation provides an approximate solution of the former 

one. In this paper, a reduction of a certain class of a system of nonlinear parabolic equa- 
Q\ ; tions du/dt = Pu + sf(u) is proposed. An amplitude equation of the system is defined 

and an error estimate of solutions is given. Further, it is proved under certain assumptions 
Qh ■ that if the amplitude equation has a stable steady state, then a given equation has a stable 

\ periodic solution . In particular, near the periodic solution, the error estimate of solutions 

^ (— | holds uniformly in t > 0. 

Keywords: amplitude equation; renormalization group method; reaction diffusion equa- 
£h tion 

: 1 Introduction 

^ ■ A reduction of a certain class of nonlinear parabolic partial differential equations (PDEs) 

o : 

-^-= e Pu + sf{u), u = u(t,x) 6 C m ,(t,x) e RxR d , (1.1) 

O ■ at 
m ; 

is considered, where s > is a small parameter, V is an elliptic differential operator with 
constant coefficient and / is a function on C m satisfying suitable assumptions. Our study 
^ ; is motivated by the following three problems. 

^ ■ Case 1. It is well known that the Swift-Hohenberg equation 

Ji = -(d 2 x + k 2 fu + s(u - u 3 ), u,xeR, (1.2) 

with a parameter k e R, can be reduced to the Ginzburg-Landau equation (amplitude 
equation) 

dA 

— = Ak 2 d 2 x A + s(A - 3A\A\ 2 ), (1.3) 
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by means of the multiscaling method @| or the renormalization group (RG) method 0]. 
Let vq be a function in some function space and give initial conditions 



w(0, x) = v ( s/sx)e ikx + v ( y[sx)e ikx , A(0, x) = v ( yfex), 

for Eqs. dl .21 ) and (11.31) . respectively. In [O, it is proved that there exists a positive number 
C such that solutions of the two initial value problems satisfy 



\\u(t, x) - (A(t, x)e ikx + A(t, x)e~ ikx )\\ < C Ve, 



(1.4) 



up to the time scale t ~ 0(1 /a) with a certain norm. In this case, a fourth-order PDE is 
reduced to a second-order PDE. 

Case 2. Let Zl = R x (0, /) be the strip region on R 2 . Consider the boundary value 
problem of a system of reaction diffusion equations on Q 



' du 




dt 




dv 




dt 




du 






y= 



(1.5) 



dv 



o, 



where /, J and k are positive constants. This type problem was introduced by Chen, Ei 
and Lin [7] to investigate a stripe pattern observed in the skin of angelfish. Under certain 
assumptions on parameters so that the system undergoes Turing instability, they formally 
derived an amplitude equation of the form 



dA 
~dt 



2d 2 



d 4 A 



+ 



(k + d)(l - d) dx 4 l-d 



(A - 3A\A\ 2 ), 



(1.6) 



without any mathematical justification. It is remarkable that the amplitude equation is a 
fourth-order equation while a given system is a second-order equation because of a certain 
degeneracy of the dispersion relation. On the other hand, a system of equations becomes 
a single equation and the number of space variables is reduced. 

Case 3. Let us consider a system of reaction diffusion equations on R 



1 du 2 ^ 

— = Dd x u + v + s(u - u ), 

dv 7 

- = Dd- x v-u, 



(1.7) 



where D > is a constant. This system can be reduced to the Ginzburg-Landau equation 



dA 9 e 
m= DdxA + 2 (A 



3A\A\ 2 ). 



(1.8) 



In this case, the order of differential equations are the same, while a system is reduced to 
a single equation. 
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A purpose in this paper is to give a unified theory of such a reduction of PDEs (11.11) . 
and give an error estimate of solutions of amplitude equations. Furthermore, we will 
partially prove a conjecture by Oono and Shiwa [GQ, which states that if a given PDE is 
structurally stable, its amplitude equation provides the qualitative features of the given 
system. For example, a stable invariant manifold of an amplitude equation implies the 
existence of a stable invariant manifold of a given system. For example, we will prove 
that Eq. dl .61) actually provides an approximate solution of the system (11.51) . Further, a 
conjecture by [|2) is solved in the following sense; if the amplitude equation (11.61) has a 
stable steady state, then the system (11.51) has a corresponding stable periodic solution. 

Since general results for (II .1) is rather complicated, we divide main results into several 
steps as follows: 

(1) In this Introduction, our main results are stated for one-dimensional problems 
u 6 C and x e R for simplicity. 

(2) In Sec. 3, the asymptotic behavior of linear semigroups generated by elliptic differ- 
ential operators are investigated. 

(2-i) Sec. 3.1 deals with the case u e C and x e R rf . The asymptotic behavior of a 
semigroup e pt is given under the assumptions (Bl) to (B3) (Propositions 3.1, 3.2 and 3.3). 

(2-ii) In Sec. 3. 2, the case u e C" and x e R' 1 is considered. The asymptotic behavior 
of a semigroup e^' is given under the assumptions (CO) to (C3) (Proposition 3.6). 

(3) Sec. 4 is devoted to nonlinear estimates and our main results are given. 

(3-i) In Sec.4.1, the case weC and x e R d is considered, which includes Case 1 above 
as an example. The definition of an amplitude equation (reductive equation) is given. An 
error estimate of solutions (Thm.4.2) is proved under the assumptions (Dl) to (D3), and 
the existence of stable periodic solutions (Thm.4.3) is proved under the assumptions (Dl) 
to (D4). 

(3-ii) In Sec.4.2, the case u e C m and x e R d is considered, which includes Case 2 and 
Case 3 above. An error estimate of solutions (Thm.4.13) is proved under the assumptions 
(EO) to (E3), and the existence of stable periodic solutions (Thm.4. 14) is proved under the 
assumptions (EO) to (E4). Theorems 4.13 and 4.14 include all previous results. 

Since we need several integers to state our final results in Sec.4, we summarize some 
of them for the convenience of the reader. 

• An integer m denotes the dimension of unknown function: u e C m . 

• An integer d denotes the dimension of space variables: x e R f '. 

• An integer M denotes the degeneracy of the dispersion relation, which determines 
the order of differentiation of the amplitude equation. For Case 1 and 3, M = 2, 
while M = 4 for Case 2. 

• An integer D (1 < D < d) denotes a dimension of the critical direction (see Sec. 3.1 
for the detail), which gives the number of space variables included in the amplitude 
equation. For Case 2, d = 2 and D = 1 . 
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• An integer N gives the number of critical wave numbers, at which the spectrum 
of the operator P is tangent to the imaginary axis. In other words, it gives the 
dimension of a center subspace (see Sec. 4.1 for the detail). An amplitude equation 
becomes a system of Af-equations. 

Although our purpose is a system of PDEs including a degenerate case as above, it 
would be better to start with a one-dimensional case for the sake of simplicity. In this 
Introduction, we suppose u e C and x e R are one-dimensional variables to state our 
main results as simple as possible. Higher dimensional problems will be treated after 
Sec.2. 

Let P(x) = ^i = Q a a^ a be a polynomial of x e R and P := P(d x ) a corresponding 
differential operator on R. For the operator P, we suppose the following: 

(Al) Re[P(i&] < for any £ e R. 

(A2) There exist w,ieR k) ± (0, 0)) and an integer M such that 

P(+ik) = +ioj, (1.9) 
P'(+ik) = ■■■ = P {M ' l \±ik) = 0, (1.10) 
P (M \ik) = P {M \-ik) * 0. (1.11) 

(A3) a q F < and P {M \ik)i M < 0. 

The assumption (Al) implies that the spectrum cr(P) of P calculated in a suitable 
space, which coincides with P(zR), is included in the closed left half plane. If cr(P) were 
included in the open left half plane, u = is linearly stable. Since we are interested in 
a bifurcation occurred at e = 0, we supposed in Eq. (| 1 .91 ) that cr(P) includes points +ia> 
on the imaginary axis. The integer M represents the degeneracy of the dispersion relation 
A = P(i£). Define the operator Q to be 

Q -?mm*L an) 

The assumption (A3) assures that P and <3 are elliptic. In this section, we further suppose 
that integers j satisfying P(ijk) = ijco are only +1 (this will be removed after Sec.2). For 
the Swift-Hohenberg equation, a» = 0, M = 2 and k is the k in Eq. dl .21) . 

Let B r : = BC' (R; C) be a vector space of complex- valued functions /onR such that 
f(x), f'(x), • • • , f (r \x) are bounded uniformly continuous. This is a Banach space with 
the norm defined by ||/|| = sup{|/(x)|, • • • , |/ (r) (x)|}. For a function / : B r — > B r , define 
the function R : (B r x B' ) -> B r to be 



R(A U A 2 ) = 



— f{A x e ikx +A 2 e~ ikx )e~ ikx dx (when k * 0), 
2tt Jo 

— f{A x e i0Jt +A 2 e-'" t )e-<" t dt (when to # 0). 



I 2n Jo 



(1.13) 



One can verify that these two expressions coincide with one another if k + and a> ± 0. 
For example if f(u) = u - u 3 , then R(Ai, A 2 ) = A\ - 3A^A 2 . 
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Now we consider two initial value problems 
du 



— =Pu + sf(u), u(0, x) = Vi(T]xy k - x + v 2 (r]x)e~ ikx (1.14) 
dt 



and 



— =QA l +sR(A l ,A 2 ), Ai(0,x) = v^x), 
dt 

2 =QA 2 + sR(A 2 ,A 1 ), A 2 (0, x) = v 2 {rjx), 



(1.15) 



dt 

where r\ : = s l/M . We call the latter system the amplitude equation. When V = -{d 2 + k 2 ) 2 
and f(u) = u - u 3 , the Ginzburg-Landau equation (11.31 ) is obtained as a special case 
A 2 - ~A\. 

Theorem 1.1. Suppose / : BC' (R; C) BC r (R; C) (r > 1) is C 1 and e > is sufficiently 
small. For any vi, v 2 e BC(R; C), there exist positive numbers C, T and to such that mild 
solutions of the two initial value problems satisfy 

\\u(t, x) - (Ai(t, x)e ikx+iajt + A 2 (t, x)e~ ikx ' ilo, )\\ < Cr] = Cs l/M , (1.16) 

for t < t < T /e. 

If we suppose A\ = A 2 := A, the system (11.151) is reduced to a single equation dA/dt = 
QA + sR(A, A); the set {A { = A 2 ] is an invariant set of (I1.15I ). Thus we put 5(A) = R(A,A) 
and consider two initial value problems 

d 

= Pu + sf(u), u(0, x) = v Y (r]x)e ikx + vi(nx)e~ ikx (1.17) 

dt 

and 

dA 

— =QA + sS(A), A(0,x) = vi(t}x). (1.18) 
dt 

For example if f{u) = u - u 3 , then S(A) = A - 3A 3 . For the equation (11.171) . we further 
suppose that 



(A4) P(iO = P(-if) and f(u) = f(u). 

That is, P{x) and f(u) are real- valued for x, u e R. In particular, if V\{nx) e R so that 
u(0, x) is real- valued, then a solution u(t, x) is also real- valued. In the next theorem, 
BC'(R; R) denotes the set of rea/-valued functions /onR such that f(x), f'(x), • • • , / (r) (x) 
are bounded uniformly continuous. 

Theorem 1.2. Suppose / : 5C r (R;R) -> BC r (R;R)(r > 1) is C 2 such that the second 
derivative is locally Lipschitz continuous. Suppose that there exists a constant <p e R such 
that 5(0) = and S'((p) < (that is, A(t, x) = (p is an asymptotically stable steady state of 
Eq. (ll.l8l) ). If e > is sufficiently small, Eq. dl . 171) has a solution of the form 

u p (t, x) = {(p + 7]ip(t, x, 77)) • 2 cos(£x + tot) . (1-19) 
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The functions if/ and u p are bounded as rj — » and satisfy 

2^/fc-periodic in x (when k £ 0) 
• 27z7<y-periodic in f (when to ^ 0), 
constant in t, x (when to = 0, k = 0, respectively), 

This w p is stable in the following sense: there is a neighborhood U c 5C(R; R) of 4> in 
5C'(R; R) such that if Vi e £/, then a mild solution u of the initial value problem (11.171 ) 
satisfies \\u(t, •) - u p (t, -)|| — > as t — > oo. 

The assumption (A4) and a space 5C'(R;R) means that this theorem holds when 
every data are real numbers. The periodic solution u p is not asymptotically stable toward 
a complex direction. These theorems are obtained as special cases of Theorems 4.2 and 
4.3 proved in Sec. 4. 

Example 1.3. For the Swift-Hohenberg equation, the estimate (11.41) immediately follows 
from Thm.1.1 by putting A 2 = A\ and v 2 = Vi. Note that the assumption for the initial 
value Vi is more relaxed than that given in [5] because we use a mild solution. To prove the 
existence of a spatially periodic solution, note that the function S (A) is given as S (A) = 
A - 3A 3 , so that A = <p = 1/ V3 satisfies the assumptions for Thm.1.2. Then, it turns out 



that Eq. dl .21) has a stable solution of the form 

2 

u p (t, x) = — cos(kx) + 0(?]), (1-20) 

which can be obtained directly without using the amplitude equation 01. 

The above theorems will be extended to more higher dimensional problems in Sec.4. 
Example 1.4. Consider the boundary value problem (|1.5t with constants /, d and k. Let 

(d 2 x u + d 2 y u) + ku-v\ ... 
d 2 x v + d 2 y v + u-v ) {LZl) 

be the linear operator which defines the unperturbed part. We suppose that d and k satisfy 
(k + d) 2 = 4d. Then, the spectrum cr(L) is the negative real axis and the origin (See Sec. 2), 
so that the system (11.51) undergoes the Turing instability when s = 0. The eigenfunction 
for 0-eigenvalue satisfying the boundary condition is given by 











■( 









<J rf )/2)M™r' c --i k -w-'--"c a22) 

We will show that the corresponding amplitude equation is given by 
(9/4 2d 2 B A A f 

T = - ik + m - d) a7 + —/ A - 3m2) - = CM) 

where rj = s 1/A . Let us consider a solution of (11.51) with the initial condition 

f =A(0,/, l ^,Ae icy +M0^x)l n * V ICV . (1.24) 
v(0,x,y) J \(k + d)/2 J \(k + d)/2 
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From theorems shown in Sec.4, it turns out that u is approximately given by 



u(t, x, y) = A(t, x)e icy + A(t, x)e~ icy + 0(rj), (1 .25) 

for to < t < T/s. Further, the system (11.51 ) proves to have a steady solution of the form 

lu p (t,x,y)\ 2 ( } / 1 \ + 0(j]X ( i.26) 

which is periodic in y and constant in x, t. The fourth-order amplitude equation (ll.231 l 
was also derived by in a certain formal way without error estimates of solutions. The 
results in this paper assure that (|1.231 i indeed provides approximate solutions and a steady 
state. 



For both examples, the spectra of the unperturbed linear operators are continuous 
spectra including the origin. Thus it is expected that when s becomes positive from 0, 
the spectra get across the imaginary axis and bifurcations occur. Unfortunately, there are 
no systematic ways to detect such bifurcations because spectrum on the imaginary axis 
is not discrete. The reduction proposed in this paper provides a systematic way to detect 
bifurcations; a bifurcation problem is reduced to that of the amplitude equation, although 
to investigate the amplitude equation is still difficult in general. 

In Sec. 2, we will demonstrate that how the amplitude equation is derived by means 
of the RG method. The RG method here is one of the singular perturbation methods for 
differential equations proposed by Chen, Goldenfeld and Oono [Q0|. In 0, it is proved for 
ODEs that the RG method unifies classical perturbation methods such as the multiscaling 
method, the averaging method, normal forms and so on. In particular, when the spectrum 
of an unperturbed linear part is discrete and a center manifold exists, the RG method is 
equivalent to the center manifold reduction; the amplitude equation gives the dynamics 
on the center manifold. This paper shows that the RG method and the amplitude equation 
are still valid even when a spectrum is not discrete and a center manifold does not exist. 
Even if there are no center manifolds, the amplitude equation provides the dynamics near 
the center subspace and it is useful to study bifurcations of a given PDE. In particular, 
Thm.1.2 means that a bifurcation may occur at s = 0; when co = 0, it is a bifurcation of a 
steady state and when a> ± 0, a /-periodic solution appears like as a Hopf bifurcation. 

Although results in this paper are partially obtained by many authors for specific prob- 
lems [HUSO, our proof is systematic which is applicable to a wide class of PDEs. From 
our proofs, it turns out that reductions of linear operators (reduction of a given differential 
operator V to <3) and that of nonlinearities (reduction of f(u) to S (u)) can be done inde- 
pendently. The reduction of linear operators is described in Prop. 3.1 and 3.6, which have 
the following significant meaning: the semigroup e p ' generated by V is approximated by 
its self-similar part. Note that the evolution equation u = Qu has the self-similar structure 
in the sense that it is invariant under the transformation 

(t, x) i-» (c M t, cx), ceR, (1.27) 

see (11.121) . Then, Prop. 3.1 implies that a won-self-similar part of e 9t decays to zero as 
t — » oo. In other words, if we apply the above transformation repeatedly, then a non-self- 
similar part decays to zero, while a self- similar part survives because it is invariant under 
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the transformation. This self-similar part defines the linear operator <3. Such a technique 
to obtain a self-similar structure is also known as the renormalization group method in 
statistical mechanics. 



2 The renormalization group method 

In this section, we demonstrate how the amplitude equation for (11.11) is obtained by the 
RG method with examples. The RG method is a formal way to find amplitude equations 
and the results in this section will not be used in later sections. Although we only consider 
parabolic-type PDEs, the RG method is applicable to a more large class of PDEs and it 
has advantages over the multiscaling method [1]. See [3] for the RG method for ODEs. 

Let us consider the Swift-Hohenberg equation (11.21 ). We expand a solution as u = 
u + eu\ + 0(s 2 ). The zero-th order term w satisfies the linear equation 

^ = -(dl + k>)W (2.1) 

We are interested in the dynamics near the center subspace. The spectrum of -{d 2 + k 2 ) 2 
intersects with the imaginary axis at the origin, and the corresponding eigenfunctions are 
e' kx and e~ lkx . Thus we consider the solution of the form 

u Q (x) = Ae ikx + Be~ ikx , (2.2) 

where A, B e C are constants. Then, the first order term u\ satisfies the inhomogeneous 
linear equation 

— = -(d x + rr«! + M - «o 

= -{d 2 x + k 2 ) 2 Ul + Ae ikx + Be~ ikx - (A 3 e 3ikx + 3A 2 Be ikx + 3AB 2 e' ikx + B 3 e~ 3ikx ). (2.3) 

Since the factors e ±lkx in the inhomogeneous terms are eigenfunctions of the operator 
-(d 2 + k 2 ) 2 , it is expected that a solution of this equation includes secular terms which 
diverge in t and x. To find secular terms arising from the factor e lkx , we consider the 
equation 

^ = -id\ + k 2 ) 2 u { + (A - 3A 2 B)e ikx (2.4) 

instead of Eq. fl2.3l) . We assume a special solution of the form 

Ml = +/ u 7 y 2 )e ikx . 



Substituting this into Eq. fl2.4l) . we obtain By = 1, B 2 = 2, and fii,/J. 2 prove to satisfy the 
relation 

ix x = Sk 2 n 2 +A- 3A 2 B. (2.5) 
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Secular terms corresponding to the factor e lkx are obtained in the same way. Therefore, a 
special solution of Eq. (l2.3l) including secular terms are given by 



B } 



A 3 

u x = (jiit + ix 2 x 2 )e lkx + (fi.it + fi 2 x 2 )e~ lkx - —e 3lkx - —e~ 3 

64 64 



where fii and ju 2 satisfy the same relation as (12.51) . Thus we obtain 

u = Ae ikx + Be~ ikx 

(nonsecular)) + 0(s 2 ). 



+ £ 



((piit + n 2 x 2 )e ikx + (fiit + p, 2 x 2 )e >kx + 



In what follows, we omit to write down nonsecular terms and 0(e 2 )-terms which will not 
be used later. To remove the secular terms, we introduce dummy parameters r and X, and 
rewrite the above u as 



u = Ae ikx + Be~ ikx + (s^r + s^ 2 X 2 )e ikx + (sfliT + sli 2 X 2 )e' ikx 
+ s[n Y (t - t) + n 2 (x 2 - X 2 )) e ikx + s(fii(t - t) + ju 2 (x 2 - X 2 )) e 



ikx 



Now terms e/^r + s/u 2 X 2 and efiyT + ejl 2 X 2 are renormalized into the constants A and B, 
respectively. Thus we rewrite u as 



■ikx 



u = A(T,X)e ikx + B(T,X)e' ikx 

+ e (jui(? - t) + ^ 2 (x 2 - X 2 )) e ikx + e (p.i(t - r) + jl 2 (x 2 - X 2 )) e 

Putting t = t and X = x provides 

u(t, x) = A(t, x)e ikx + B(t, x)e~ ikx , 

which seems to give an approximate solution if A(t, x) and B(t, x) are appropriately de- 
fined. Since u is independent of dummy parameters t and X, we require that the equation 



du d 2 u 







3t dx 2 

holds, which is called the RG equation. This yields 



du 
d 2 u 



T=t,X=X 



OA 



OX 2 



T=t,X=X 



dt 
d 2 A 
dx 2 



= \^--Efii\ e ikx + \^--efii\ e~ lkx = 0, 



dB 



dt 



2six 2 \e ikx + 



d 2 B 

dx 2 



(2.6) 



2sfi 2 \e- ikx = 0. 



Since [i\ and fx 2 satisfy (12.51) . we obtain 
OA 



d 2 A 

^-=s/ui= Sk 2 su 2 + s(A - 3A 2 B) = 4k 2 — 
dt ox 1 



+ s(A - 3A 2 B). 



(2.7) 



Similarly, B satisfies dB/dt = Ak 2 d 2 x B + s(B - 3AB 2 ). If we suppose that A = B to obtain 
a real-valued solution, the Ginzburg-Landau equation (11.31) is obtained. 
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Next, let us derive the amplitude equation of the system (11.51) . The dispersion relation 
of the unperturbed operator L (11.211) is 



det 



-1 A + + 1 



A 2 + (df + £ 2 + 1 - k)A + (d£ 2 - k)(£ 2 + 1) + 1 = 0, 



(2.8) 



where we put £ 2 = ^ 2 + Let A±(^) be two roots of (12.81 ) . Then, the spectrum of L is 
given by <x(L) = A+(R) U /L(R). Suppose that Eq. d 1 .51) undergoes the Turing instability at 
s = 0, so that cr(L) = R< . It is easy to verify that this is true if and only if 

< d< k < 1, (k + d) 2 = Ad. 

In particular, one of A ± (g) satisfies A+(c) = 0, where c 2 = (k - d)/2d. For any (£1,^2) 
satisfying = c 2 , 



1 



(k + d)/2 t 

is an eigenfunction of L associated with /I = 0. Because of the boundary condition in 



131) . we choose e lcy and e '° . Thus we expand a solution of (11.51) as 



1 

(Jfc + d)/2 



+ B 



1 

(it + d)/2 



e" cy + e 



+ 0(s 2 ). 



(2.9) 



Put A = B for simplicity. Then, (u\, V\) satisfies the equation 







!::) 




(::) 


-L 


+ ( 











Ae icy _ A 3 e 3icy _ 3A \ A \2 e icy + c ( . 





where c.c. denotes the complex conjugate. We find secular terms of the form 

k + d 

u\ = (jj.it + H2X 2 + /j.3X 4 )e' cy , vi = — - — (fiit + fax 2 + fax A )e' cy . 



(2.10) 



(2.11) 



Substituting them into (12.101) . we obtain 



A*i = Mi, 

fa = i"3, 
L fa = 2^1 > 



fii = 2dfj. 2 + A- 3A\A\ 2 , 

k + d 2 
= 12J// 3 + kfi 2 r — fa - dc /J.2- 



(2.12) 



Then, a formal solution is given as 

u = Ae lcy + s(juit + H2X 2 + /J3X 4 )e' cy + c.c. + (nonsecular) + 0(s 2 ). 
Introducing dummy parameters r, X and renormalizing, we rewrite this equation as 
u=A(T,X)e icy +s(m(t-T) + ju 2 (x 2 -X 2 )+iu 3 (x 4 -X A ))e ky +c.c.+(nomecu\ar) + 0(£ 2 ). 
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Since u is independent of r and X, we require 



( du 

d 4 u 



T=t,X=X 



dX 4 



OA 
d 4 A 



— -sy. x \e } +c.c. 



0, 



T=t,X=X 



dx 4 



(2.13) 



-24eiu 3 \e icy + c.c. = 0. 



Finally, Eqs. (l2. 121) and (12.131) provide the amplitude equation (11.231) by eliminating iu.i,/u 2 
and yu 3 . 



3 Reduction of a linear semigroup 

For Eq. Ol.ll) . reductions of the linear unperturbed part Vu and the perturbation term f(u) 
can be done independently. In this section, we give a reduction of the linear part. 



3.1 One dimensional case 

We start with the simple case u e C and x e R f '. Put x = (xi, • • • ,x d ) and a = (ai , ■ ■ ■ , a d ), 
where a denotes a multi-index as usual: x a = (x" 1 , • • • , x a d d ) and \a\ = ot\ + ■ ■ ■ + ad- Let 
P(x) = Y!l a \ = Q^aX a be a polynomial of degree q and V := P{d\, ■ ■ ■ ,dj) a differential 
operator on R J , where dj denotes the derivative with respect to Xj. We make the following 
assumptions. 

(Bl) Re[P(i£)] < for any £ e R d . 

(B2) There exist co e R, k e R f/ and an integer M such that 
P(ik) = iu), 

gap 

— — (ik) = 0, for any a such that \a\ = 1, • • • , M - 1, 

ox" 

d a P 

(ik) ± 0, for some a such that lal = M. 

dx" 

(B3) Define Q(x) and (3 by 

= X r-irrrSw^ Q = ■ ■ ■ • ^- C 3 - 1 ) 

Then, both of !P and (3 are elliptic in the sense that there exist C\,C2 > such that 
Re[P(i#] < -c 2 |£| 2 and Re[G(i#] < -c 2 |£| 2 hold for |f| > Cl . 

Put 5' = 5C r (R rf ; C), a Banach space of complex-valued bounded uniformly continu- 
ous functions on R'' up to the r-th derivative. In the next propositions, || • || = || • || ( . denotes 
the standard supremum norm on B r . Consider two initial value problems: 

u(0,x) = v (x)e ikx , (3.2a) 



A(0,jc) = v (jc), (3.2b) 
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where kx = k\X\ + ■■■ + k^Xd, and a similar notation will be used in the sequel. Because 
of (B3), V and <3 generate C°-semigroups e 91 and e Qt on B r , respectively. Thus solutions 
of the above problems are written as e p '(e' kx v ) and e Qt v , respectively. 

Proposition 3.1. Suppose (Bl) to (B3) and r > 0. There exists a constant C\ > such 
that the inequality 

||e*V*v ) - e Mx e Qt v Q \\ r < C { r XIM \\v Q \\ r (3.3) 

holds for any t > and v 6 BC r (R d ; C). 

For the main theorems in this paper, we need the following perturbative problem 

( du 

— = Pu, u(0,x) = v (r]x)e lkx , (3.4a) 
at 

dA 

— =QA, A(0,x) = Vo(t]x), (3.4b) 
I at 

where r\ = s i/M and e > is a small parameter. 

Proposition 3.2. Suppose (Bl) to (B3) and r > 1. For any e > and ? > 0, there exists a 
positive number Ci = C\(to) such that the inequality 

\\^(e ikx v Q ) - e^ +ikx e Q %\\ r < r,C { \\v Q \\ r (3.5) 

holds for t > t and v 6 BC(R d ; C), where v (0 := v (v ■)■ 

Proof of Prop.3.1. By putting m = e'^w, Eq. d3.2at is rewritten as dw/dt = (P - iio)w. 
Then, the operator V - icj satisfies (Bl) to (B3) with a> = 0. Hence, it is sufficient to prove 
the proposition for co = 0. 
Two solutions are given by 



and 



u(t,x) = -JL Jv (y + x)e ik( y +x) Je^e^'d^dy 

= (^/ v o^ + ^) fe-Me^^dtdy, (3.6) 
respectively. Thus we obtain 

u(t,x) - e ikx A(t,x) = J v Q (y + x) J e Qm (e P{i ^ ik)t - Qm - l)d£dy. 
Put r = r 1/M . Changing variables £ i-» t£, y i-» y/r yields 

x) - e fa A(f, x) = J v (y/r + x) J - l) dfdy. 
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Due to the assumption (B2), we have 



g(£,T) := P(iT{ + ik)/T M -Q(i{) 



= -y 

M=o 

z 



i 



(ar 1 !)---(or rf !)a* 



\a\=M 



( ai \)---(a d \) dx« 



1 



gap 



(anl)- ■■(a d \)dx a 



\a\ga -\a\-M 



\a\=M+l 

Note that g ~ 0(t) as t — > 0. In particular, there exists < 6 < 1 such that 



^ r) -l=T|(f,9r) 

OT 



This provides 



n(f, x) - e ikx A(t, x) = t 



(2 



v (y/T + x)G(y, r)dy, 



G(y,r):= f e-^e Q(i ^(Z,6T)e^ 6T) dt. 
I J dr 



(3.7) 



Because of (B3), G(y, r) exists for each r > and v e R. Since g is polynomial in 
r, there exist To and D\ = Di(to) such that |G(y, r)| < D\ holds for < r < To and 
y e [-1, l] d . Next, since the integrand in the definition of G(y, r) is smooth in G(y, r) 
is rapidly decreasing in y due to the property of the Fourier transform. Indeed, by using 
integration by parts, it is easy to verify that there exists D 2 = D 2 (tq) such that \G(y, t)| < 
Z) 2 Cyi • • - yd) 2 holds for < r < To and y [-1, l] rf . This provides 



«0, jc) - e ikx A(t, x)\ < 



(2 



|v (y/T + x)|-|G(y,T)|Jy 



< 



llvoll + 
(2n) d (2. 



r_ r _jh 

*) ci Jyt[-i,u*y 2 ■ ■ ■ 



dy 



This proves that 



sup \u(t, x) - e ikx A(t, x)\ < tD 3 | |v | | 

xeR d 



(3.8) 



for some D 3 > when < r < r . To estimate the derivatives, note that Eq. fl3.6l) is 
rewritten as 



e ikx u(t,x) 



(2 



l — d j v (y) j e^e^d^y, 



and similarly for A(t, x). Hence, the derivative is given as 

( d a 



— (f ikx u(t, x) - A(t, x)j = t—^ J v (y/r + x)G a (y, r)dy, 
G a (y, r) := J 'jT&e-Ke*®^, 6r)e^d^. 
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By the same way as above, we can show that this derivative is of 0(r) uniformly in x. 
Hence, the inequality 

\\u(t,x) - e ikx A(t,x)\\ < tD 3 \\v \\ = r 1/M D 3 ||v || (3.9) 

holds with respect to the norm of B' for some D 3 > and any t > t^ m . On the other hand, 
since V and <3 generate C°-semigroups on B r , there exists D 4 > such that \\u(t,x) - 
e ikx A(t, x)\\ < D 4 \\v \\ for < t < t q m . This and Eq.(TJU) prove Prop.3. 1 (for to = 0). ■ 

Proof of Prop.3.2. In this case, solutions satisfy 

g ikx r 

u(t, x) - e lkx A{t, x) = t— — voirjy/T + rjx)G(y, r)dy, 

(2n) d J 

where t = r i/M and G is defined by (13.71) as before. Since v e B' (r > 1), there exists 
< Q\ < 1 such that it is expanded as 



£ ikx 

u(t, x) - e' kx A(t, x) = t 



r( d dv 

I v (t]x) + V — -(Tjx + etTjyM-yj 

J [ OXj T 



(2n) d 

I V —^-(nx + 



G(y, T)dy 



11 '{2nyi J lj-fa(Vx + OiriyMyjG(y,T)dy, 

7=1 3 



where we used the fact fG(y, r)dy = 0. The rest of the proof is the same as that of 
Prop.3.1. ■ 

If the polynomial P(x) has no symmetries, the assumptions (B2),(B3) seem to be 
strong; for example, if d = 2 and 



d 2 P d 2 P d 2 P 

—m * o, - — (ik) = — 

ax\ oxix 2 ax\ 



then Q is not elliptic. To relax the assumptions, fix an integer D such that 1 < D < d. We 
denote x e R' 1 as x = (x\, x 2 ) with x,\ = (xi, • • • , x D ) and x 2 = (x D+ i, • • • , x t i). Accordingly, 
a multi-index a is also denoted as a = (J3, y). Instead of (B2) and (B3), we suppose that 

(B2)£» there exist co eR, k eR d and an integer M such that 

P(ik) = ia), 
&P 

— -(ik) = 0, for any B such that \B\ = 1, • • • , M - 1, 
dx { 

d?P 

— -(ik) ± 0, for some B such that \B\ = M. 



(B3) D Define Q(x) and Q by 



Q(x) = Q(x x , 0) = ^ — — — — — -(ik)x\, Q=Q(d u ---,d D ,0,---,0). (3.10) 
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Then, both of V and <3 are elliptic in the sense that there exist Ci,c 2 > such that 
Re[i'(^)]<-c 2 |^| 2 andRe[e0| 1 )]<-c 2 || 1 | 2 holdfor|^|,|| 1 | > c u where ft = (£,••• ,% D ). 

When D = d, these assumptions are reduced to (B2) and (B3) before. The assumption 
(B3)d implies that Q is an elliptic operator on R D although it is not on R rf . Consider two 
initial value problems: 



( du 

— =Pu, u(0,x) = v (x 1 )e lkx , 
Ot 

dA 

— =QA. A(0,x) = v (x 1 ). 
K ot 



(3.11a) 



(3.11b) 



Note that vo depends only on ii = (x\,--- ,x D ). In particular, Eq. (l3.1 lbl) can be re- 
garded as a parabolic equation on BC r (R D , C), while Eq. (l3.1 lal) is a parabolic equation 
on BC'(R d , C). We also consider the perturbative problem 



■^=Pu, u(0,x) = v (7 1 x l )e ikx , 
ot 

dA 

— =QA, A(0,x) = v (t]x 1 ), 
y ot 



(3.12a) 
(3.12b) 



where r\ = s i/M and e > is a small parameter. 

Proposition 3.3. Under the assumptions (Bl), (B2) D and (B3)£>, Prop. 3.1 holds for 
Eqs. (l3TTTaT) . (l3TTTbl) . and Prop.3.2 holds for Eqs. (l3TT2^ . (l3T2bl) . 

Proof. For Eq. d3.Hal l. u(t, x) is given as 



u(t,x) = Jvo^i+xO Je-^e p ^ +ik)t d^dy 



e 

(2 
e 
(2 



ikx p n 



iyi f i -iyi& P(i%\ +iki ,ifr J rik 1 



^d^id^jdyidh 



To calculate this, we need the next lemma. 

Lemma 3.4. Let S be a space of C°° rapidly decreasing functions on ~R d ~ D (Schwartz 
space). For any / e S, we have 



m)dhd% = (2n) d - D f(0). 



(3.13) 



Proof. Let S' be a dual space of S. For the pairing of S' and <S, we use a bracket ( , ). 
Let T be the Fourier transform. Then, 



(d-D)/2 



r[f](y 2 )dy 2 



= (2nf- D)/2 (l,T[f]) 
= (2ny d - D)/2 (T[\],f) 
= (2n) d ' D (6,f) = (2n) d - D f(0), 
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where 6 is the Dirac delta. 



Due to this lemma, we obtain 



Since Eq. (l3.11bl) is an equation on BC r (R D , C), A(t, x) is given as 
The rest of the proof is the same as those of Prop. 3.1 and 3.2. 



3.2 Higher dimensional case 

Suppose w = («i,---, u m ) 6 C m and x = (xi, ■ • • , Xd) e R d . For fixed 1 < Z) < d, we use 
the same notation x = (jci, x 2 ) as in Sec.3.1. Let {Pij(x)}f ]=l be the set of polynomials of 
x. Define the matrix P(x) by 



P(x) = P(x u - • • ,x d ) 



Pn(x) 



Plm(x) 



Pml(x) 

The differential operator V is defined to be V = P(d\, 

(A-P n (iO ■■ 



(3.14) 



det(/l - P(g)) = det 



Pmm\X) , 

, <9d). The algebraic equation 

-Pi m m ^ 



A ~ Pmm(ig) 



= 



(3.15) 



is called the dispersion relation. Let /li(^), • • • ,A m (g) be roots of this equation. Then, 
/ti(R) U • • • U A m (R) gives the spectrum of V. We suppose for simplicity that only A\{£) 
contributes to the center subspace of V (see (CI) below). Extending to more general 
situations is not difficult (see Remark 3.7 below). 

(CO) The matrix P(i£j) is diagonalizable for any ^ e R''. 

(CI) Re[Ai(#] < and Re[Aj(£)] < for any f 6 R d and j = 2, • • • , m. 

(C2) There exist weR^eR" 1 and an integer M such that 

A\{k) = ico, 
eft A 

-(k) = 0, for any fi such that = 1, • • • , M - 1, 



as? 



(fc) # 0, for some fi such that \J3\ = M. 
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(C3) Define Q(x) and Q by 
Q(x) = Q(x u O) = Yj 



= M 



(ft !)•••<&,!) dt 



(k)( Xl /if, Q = Q(d u - - ,d D ,0,--- ,0). 



(3.16) 

Then, both of f and <3 are elliptic in the sense that there exist c\,ci_ > such that 
Re[A/#] < -c 2 |^| 2 (7 = I,--- ,m)andRe[Q(/| 1 )] < -c 2 |#! | 2 hold for |£|, &I > Cl , where 
|i = (fi, •••,&). 



Put 5 r = BC'(R r ; C) and let (B r ) m be the product space. The norm on (B r ) m is defined 
by ||w|| = max^j^,, \\uj\\ for u = (wi, • • • , w m ). Note that "P is an operator densely defined 
on (B r ) m while (3 is an operator densely defined on B r . When m = 1, A\(£) = P(i£), so 
that the above assumptions and Q are reduced to those in Sec.3.1. 

Example 3.5. Suppose m = d = 2 and consider the operator L defined by (11.211) with 
the condition < d < k < 1, (k + d) 2 = Ad. The dispersion relation is given by (12.81 ). 
whose roots are denoted as A 2 {^) < Ai(£). It is easy to verify that A\(g) = if and only if 
f = (£1,^2) satisfies £ 2 + %\ = c 2 := (k - d)/2d. Thus there are infinitely many points £ 
satisfying A\(g) = 0. We choose (£1,^2) = (0, c). Then, we can show that 

dX\ ,„ . d 2 A x , n x d^Ax , n x rt d 4 ^ 



Ai(0,c) = — -(0,c) 



2 (0, c) = — ^(0, c) = 0, — ^(0, c) # 0, 



while 



|i(0,c) = 0, ^(0, 



Hence, (C2) is satisfied with k = (0, c), co = 0, D = 1 (i.e. and x 2 = y), and M = 4. 

In this case, 

2d 2 <9 4 



(3 



1 3% ,„ . d 



(0,c)— = - 



(3.17) 



4! dtf^'^dx 4 (k + d)(\ - d) dx*' 
see Eq.CC23]>. 

Let w = (wi, • • • , w m ) be an eigenvector of P(ik) associated with A\(k) = ito. Note that 
e lkx w is an eigenfunction of V included in the center subspace. Consider two systems of 
PDEs: 



- r =Pu, u(0,x) = v (x l )e ikx w, 
Ot 

■4= OA. A(0,x)=v (xi). 



(3.18a) 



(3.18b) 



Note that v depends only on x\ = (x\, ■ ■ • , x D ). The former is a system of /^-equations on 
R (/ , while the latter is a single equation on R D . We also consider the perturbative problem 



du ■, 

— = Pu, u(0, x) = v (rjxi)e l x w, 
ot 

dA 

— =QA, A(0,*) = v (J7*i), 
I ot 



(3.19a) 
(3.19b) 
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where 77 = s 1/M and s > is a small parameter. Solutions of them satisfy the next 
proposition. 

Proposition 3.6. Suppose r > 0. Under the assumptions (CO) to (C3), there exists a 
constant C\ > such that 

||/V* v vo • w) - e' w+,fa (e Q V ) • w|| < C { r llM \\vo\\ (3.20) 

holds for any t > and v e 5C'(R' ; C). Next, suppose r > 1. For any £ > and ? > 0, 
there exists a positive number Ci = C\(t ) such that the inequality 

lk*V**v • w) - *> m *Vv ) • w|l < ^Cillvoll (3.21) 
holds for t > t and v e 5C'(R r ; C), where v (-) := v (r] •). 

Proof. We suppose D = d for simplicity; that is, k\ = x and f3 = a. The case D < d 
is easily reduced to the case D = d as in the proof of Prop. 3. 3. We also suppose oj = 
without loss of generality. 

Like as the proof of Prop. 3.1, a solution of (I3.18al) is written as 

u(t, x) = ^ d jv (y + x) J e**e™***wd&y. 
Note that e p ^ +,k)t is an exponential of a matrix. Let S (£) be a matrix such that 

(HO 



S(&- l P(.i&S(.fi:=A(& 



(3.22) 



Because of the assumption (CO), we can assume that S (£), S (£) 1 and /l/£)'s are smooth 
in Then, 

Put t = r 1/M . Changing variables £ i-» t£, y i-» y/r yields 

M(f ' X)= (2^ Jv (y/r + x) Jg-^S^ + T^^^S^ + T^wdfify. 
Expanding S(£ + r£) and 5(fc + r^)" 1 , it turns out that there is a function G\ such that 

u(t, x) = ^ J i/„(y/T + x)J r^S (^^^^S^-'W^ + r J Vo (y/r + x)G l (y, r)dy. 

By a similar estimate used in the proof of Prop.3.1, we can show that there exists D\ > 
such that the norm of the second term above has an upper bound tDi||vo|| for any r > 0. 
Since w is an eigenvector associated with A { (^), we obtain 

u(t, x) = J v (y/r + x) J e^ k+r ^ wd^dy + rj ^(y/r + x)G l (y, r)dy. 
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Therefore, we have 



u(t, x) - e ikx A(t, x)w 



(2? 

where g(£, r) = A Y (k + t£)/t m - Q(i^). The rest of the proof is the same as those of 
Prop.3.1 and 3.2. ■ 

Remark 3.7. Let A\(^), ■ ■ ■ ,A m (^) be eigenvalues of P(i<;) as before. Even if several 
eigenvalues lie on the imaginary axis and (CI) is violated, to modify Prop. 3. 6 is very easy; 
since equations are linear, the superposition principle is applicable. A typical problem is 
that P(ig) is real-valued and eigenvalues occur in complex conjugate pairs. For example, 
suppose (CO), (C2), (C3) and the following (CI)' instead of (CI): 



(CI)' A 2 (0 = Aitf). Re[A l2 (£)] < and Re[Aj(g)] < for any £ e R' and j = 3, • • • , m. 



Put Q := Q(d\, •■• , d D , 0, • • • ,0). Let h>i and w 2 be eigenvectors of P(ik) associated with 
A\(k) = ia) and A 2 (k) = -ico, respectively. In this case, instead of Eq. fl3.20l) . the inequality 

||c^(«**vi • wi + e ikx v 2 ■ w 2 ) - e^^ie^Vi) ■ w x - £T iM *(e S V 2 ) -w 2 \\< C x t~ l ' M {\\v { \\ + \\v 2 \\) 

(3.23) 

holds for any t > and V\,v 2 e BC r (W; C), and similarly for Eq. (l3.21l ). 
Example 3.8. Suppose m = 2 and d = 1 . Define a linear operator 

u \ I Dd 2 1 \ u 



P U= -1 D# ,1- < 3 ' 24) 



where D > is a diffusion constant. This operator arises from Eq. (11.7t . Eigenvalues of 
P(i£) are A x (t-) = -D% 2 + i and A 2 (%) = ~D% 2 - i. Hence, the assumptions (CO), (CI)', 
(C2) and (C3) are satisfied with 

d 2 

AM = i, A 2 (0) = -i, M = 2, Q = Q = D-, 
W\ = ( 1 ), w 2 = | . | = w,. 



-I 

Eq. (l3.23l) is given as 

\\e p \ Vl -w l+ v 2 - w 2 ) - e i0Jt (e Qt Vl ) ■ w, - e^W) • w 2 \\ < dr 1/Af (||vi|| + ||v 2 ||). (3.25) 



In most applications, we take v 2 (x) = Vi(jc) to obtain a real-valued solution of ii = Vu. The 
above inequality implies that an approximate solution of u = Vu is constructed through 
the complex heat equation A = OA. 

4 Main theorems 

In this section, a reduction of a perturbation term f(u) is given. Combined with the 
reduction of linear semigroups, a reduction of Eq. dl . It is performed. 
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4.1 One dimensional case 

We start with the case u e C and x e~R d . Put x = (xi, • • • , xj) and a = {a\ , • ■ • , o^), where 
a denotes a multi-index. For a fixed integer 1 < D < <i, we denote x e R d as x = (xi, x 2 ) 
with xi = (xi, • • • , x D ) and x 2 = (x D+i , ••• , x d ). Accordingly, a multi-index a is also 
denoted as a = (J3,y). Let P{x) = ZJ^q fi^x 0, be a polynomial and !P := P(d\, ■ ■ ■ ,dd) 
a differential operator on R rf , where <9 ; - denotes the derivative with respect to x y . For the 
main theorems, we make the following assumptions. 

(Dl) Re[P(i£)] < for any £ e R rf . 

(D2) There exist co e R, k e R rf ((w, fc) ^ (0, 0)), a finite set of integers / = {ji, • • • , 
and {Mi, • • • , M N ) such that 

^07n&) = ijnio, (n= 1, • • • , A 7 ), 

— -(ij n k) = 0, for any 8 such that |/?| = 1, • • • , M„ - 1, (n = 1, • • • , AO, 

— oiijnk) ± 0, for some 8 n such that \8 n \ = M„, (n = 1, • • • , AO- 

dx? 

The set J consists of all integers satisfying P(ijk) = ijco. 
(D3) For n = 1, • • • , N, define Q n {x) and Q n by 

Zl c^P 
73-77 77777 — j(ijnk)f v Qn = Q n {du ' ' - , d D , 0, • • • , 0). 

(4.1) 

Then, both of V and <3„ are elliptic in the sense that there exist c\,C2 > such that 
Re[P(m < -c 2 \£\ 2 and Re[Q„(i|i,0)] < -c 2 \U hold for |£|, &I > c, and n = 1, • • • ,N, 
where |i = (ft,-- - ,£ D ). 

In addition to (B2) before, a new assumption fc) ± (0, 0) and the set J are intro- 
duced. In most examples, / consists of / = {+1, -1} as Sec.l, see also an example below. 

Put B' = 5C"(R d ; C). For a given function / : B r — » 5'", let us consider the Fourier 
series of the quantity /(X^ =1 A n e i] ' n0>t+ijnkx ), where A = (A x , • • • , A w ) e C w . Since (w, fc) 
(0, 0), the Fourier series is well-defined and it is easy to verify that the series is of the form 

N co 

fiJ^Ane'j"^"^) = ^ Cj{A)e i3 " t+ijkx . (4.2) 

n=l j=~°° 

For example, when k\ ± 0, C 7 (A) is given by 

Cj(A) = Cj(A u ■ • • ,A N ) := -i- /( V Aj^e'^ dx x . (4.3) 

When a» ^ 0, it is also written as 

C ; (A) = — /(VA/^V-^. (4.4) 
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In particular, Cj n (A) is denoted by R,,(A) if j„ e J. For any 6 e R and j e Z, Cj(A) satisfies 
the equality 

Cj(e ihe A u - ■ ■ , e i]N6 A N ) = e i]e Cj(A u ■ ■ ■ ,A N ). (4.5) 
Let s > be a small parameter. We will consider the two initial value problems: 



du 



N 



— = Vu + sf(u), u(0, x) = ^ e^v^), 



dt 
dA 



n=\ 



v dt 



- = Q n A n + eR n {A), A„(0, x) = v n (T]Xi ), (n =!,-■■ ,N), 



(4.6a) 



(4.6b) 



where 7] = s l/M and M := min{Mi, • • • ,M N }. Note that the former is a single equation 
while the latter is a system of PDEs. 

Example 4.1. Let us consider the Swift-Hohenberg equation (I1.21 ). For this equation, 
D = d = 1 and />(*) = -(x 2 + k 2 ) 2 . Since = if and only if £ = the set J consists 
of ji = + 1 , j% = - 1 . We have 

— (±ifc) = 0, -^ji+ik) = 8k 2 . 
ox our 

Thus M\ = M 2 = 2, and both of Qi and <3 2 are given by 

Q -l^L( +ik )—-Ak 2 — 

2 dx 2 dx 2 dx 2 

Since f(u) = u - u 3 , the expansion of f(Aie lkx + A 2 e~ lkx ) is 

f{A x e ikx + A 2 e~ ikx ) = A x e ikx + A 2 e~ ikx - {A\e 3ikx + 3A 2 x A 2 e ikx + 3A\A\e~ ikx + A 3 2 e~ 3ikx ). 
This provides 

tfi(A) = Ci(A) = Ai - 3A 2 A 2 , # 2 (A) = C_i(A) = A 2 - 3AiA 2 . 



Therefore, the amplitude equation (I4.6bl l is given by 



dAi -,d 2 Ai r, 
= Ak 2 ^TT + s(A l -3A 2 A 2 ), 



dA 2 i d 2 A 2 , 
_l = 4^ + E (A 2 -3A^). 



& "" dx 2 

Usually, we assume Ai = A 2 , which gives the Ginzburg-Landau equation (fT3 



(4.7) 



Put v„(x) = v n (i]x). The above equations are rewritten as integral equations of the form 



u = e p, {Y e ij " kx %) + e e nt ~ s) f(u(s))ds, 

f e Qn(ts) Rn 

Jo 



e~" v„ + s 



(A(s))ds, (zi=l,... , AO, 



(4.8a) 
(4.8b) 
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whose solutions are called mild solutions. When / : B r — » B' is C 1 , then R n : B r ^> B r 
is also C 1 , and due to the standard existence theorem (see Pazy[|8)), there exists a positive 
number T$ > such that the above integral equations have mild solutions u(t, •), A n (t, •) e 
B r for < t < Tq/s. In particular, when the initial condition {v n } is included in the domain 
of V and <3„, then a mild solution is a classical solution which is differentiable in t > 0. 
In this paper, we only consider mild solutions. The main theorems for a one-dimensional 
case are stated as follows: 

Theorem 4.2. Suppose / : BC r (R d ;C) BC r (R d ;C) (r > 1) is C 1 and s > is 
sufficiently small. For any {v„}^ =1 c 5C'(R rf ; C), there exist positive numbers C, T Q and t Q 
such that mild solutions of the two initial value problems (4.5) satisfy 

N 

\\u(t, x) - ^ A„(f, x)//"-'+''.>»^|| < c»7 = Ce 1/M , (4.9) 
n=i 

for ? < * < T Q Is. 



Next, let us show that the error estimate above holds for any t > to under a suitable 
condition. For ordinary differential equations, it is proved in [0 that if the amplitude 
equation has a stable hyperbolic invariant manifold, then a given equation has a stable 
invariant manifold of the same type and approximate solutions are valid for any t > 
near the manifold. For our situation, suppose that there is a constant vector (p £ R w such 
that R{cp) = 0, R = • • • ,Rn)- Then, is a steady state of the amplitude equation. 
Unfortunately, (f> is not hyperbolic because of the symmetry (14.51 ); the Jacobi matrix of R 
at (p has a zero-eigenvalue in general. For example, although the amplitude equation (14.71) 
for the Swift-Hohenberg equation has a steady state (A\,A2) = (1/ V3, 1/ V3), the Jacobi 
matrix of R at (1/ V3, 1/ V3) has a zero-eigenvalue. However, if we restrict solutions to 
the invariant set {Ai = A 2 }, (14.71) is reduced to 

d A = 4k 2 ^ + s(A - 3A 3 ), (4.10) 

and the derivative of the function A - 3A 3 at A = 1 / V3 is negative. This implies that 
A = 1 / V3 is a hyperbolically stable steady state of (14.101) . and we expect that the Swift- 
Hohenberg equation also has a corresponding stable solution. For more general situations, 
we make the following assumption. 

(D4) For \B\ = M n andn = 1, • • • ,N, 

P(i& = PHl), f(u) = W) and ^(Unk) = ^("^)- 

The first two equalities imply that P(x) and f(u) are real- valued when x, u e R. Due 
to this assumption, P(-ij n k) = -ij n co when P(ij n k) = ij n (x>. Hence, the set J of integers 
satisfying P(ijk) = ija) is given by J = {j u ■ ■ ■ , j N } U {-j u ■■■ , -j N }. We denote -j n by 
j- n . Then, M n and <3„ are defined for n = +1, • • • , +N as in (D2), (D3). It follows from 
(D4) that M n = M_„ and Q n = <3_„. For many examples, J consists of 7 = {+1} and 
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this assumption is satisfied. In the present notation, the function Cj(A) defined by (14.31) is 
given by 

Cj(A u --- ,A N ,A-u--- ,A- N )= + f /(^A„ e ^ Jcl +^A_ ne -^ j;i ) e -^^Jx 1 .(4.11) 

n=\ n=l 

In particular, Cj n (A) is denoted by R n (A) for n = +1, • • • , ±N. Hence, the amplitude 
equation is given as a system of 2./V-equations of the form 

dA 

= Q n A n + sRM), (n = ±1, • • • , ±N). (4.12) 

ot 

This system can be reduced as follows; It is easy to verify that Cj satisfies 

Cj(A u - ■ ■ ,A N ,A- U ■ ■ ■ ,A- N ) = C„j{A„ u - ■ ■ ,A_ N ,A U ■ ■ ■ ,A N ). (4.13) 

Thus putting A n = A_„ yields R n (A) = i?_„(A). Since Q n = Q „, putting A„ = A_„ shows 
that Eq. (l4.12l) is reduced to the system of iV-equations 

dA 

= Q„A n + sR n (Ai, ■ ■ ■ ,A N ,Ai, ■ ■ ■ ,A N ), (n = l,---,N). (4.14) 

ot 

Define the function S „ to be 

S n (A) = S n (A u --- ,A N ) = R n (A u --- ,A N ,Au--- ,A N ). (4.15) 
We consider the two initial value problems: 

du 



— = Vu + sf(u), «(0, x) = J] {e iJnkx + e~ iJnkx ) v B (^i ), (4.1 6a) 

n=l 

dA 

= Q n A n + sS„(A), A„(0, x) = Vnirjxt), (n = l,---,N), (4.16b) 

v at 

where rj = s 1/M and M := min{M 1; • • • , M N }. Due to Thm.4.2, solutions of them satisfy 

N 

\\u(t, x) - J^Anit, X ){e ik0,1+ii " kx + e-'J-w-'J-^w < Crj, 

for t < t < T /e. Further, we can show the next theorem, in which B r = 5C(R f ';R) 
denotes the set of real-valued functions. 

Theorem 4.3. Suppose (Dl) to (D4) and / : BC r (R d ;R) -> 5C(R f/ ;R) (r > 1) is C 2 
such that the second derivatives are locally Lipschitz continuous. Suppose that there exists 
a constant vector <f> = {<p\, ■ ■ ■ , <p N ) e R N such that 

(i) S n ((f>) = Oforn = l,--- ,N, 

(ii) the Jacobi matrix of (Si, - •• ,Sn) at (f> is diagonalizable and all eigenvalues of the 
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matrix have negative real parts. 

If e > is sufficiently small, Eq. (l4.16al) has a solution of the form 



,($, x) = 2_j(<f>n + rj^nit, x, 77)) • (e' Jl 



a)t+ij„kx + e -ij„ajt-ij n kx\ 



(4.17) 



The functions tff n and u„ are bounded as rj — » and satisfy 



27r/o»-periodic in t (when co ± 0), 
constant in t (when co = 0), 



27z7& ; -periodic in x 7 - (when kj 0), 
constant in x; (when kj = 0), 



for j = 1, • • • , J. This w p is stable in the following sense: For any n = 1, • • • ,7Y, there is 
a neighborhood [/„ c BC r (R d ; R) of 0„ in BC(R d ; R) such that if v n e U n , then a mild 
solution m of the initial value problem (14. 1 6ab satisfies \\u(t, •) - u p (t, -)|| — » as f — > 00. 

The above conditions (i),(ii) show that A(t, x) = cp is an asymptotically stable steady 
state of Eq. (l4.16bl) . Thus this theorem implies that a stable steady state of Eq. (l4.16bl) 
induces a periodic solution of Eq. (l4.16al) . Due to the symmetry (I4.51 ). Eq. (l4. 121) has a 
steady state A n = e lj " 6 (p n (n = +1, • • • , +A0 for any 8 e [0, 2n). Accordingly, we can prove 
that Eq. (l4.16al) has a stable periodic solution 

N 

u p (t,x) = Yjpn + Tp/f n (t,x,jf)) ■ (e ij » ( ° t+ij " kx+ij » e + e-wt-ijnkx-Hne^ 

n=\ 

for any 6 e [0, In), the proof of which is reduced to (14.171 ) by the translation of t or x. 

Proof of Thm.4.2. We prove the theorem for the case D = d (thus x.\ = x and a = B) for 
simplicity of notation. The general case D < d can be proved in the same way. A proof is 
divided into four steps. 

Step 1. notation. It is convenient to introduce some notation: We define a new coordi- 
nate (T,X) by 

x = X/tj, t = T/s, u(T,X) = u(t, x), A n (T,X) = A n (t, x), 

P = -P(r]d x ), 4 = -Qnirjdxl 
s e 

Then, Eqs. (l4.6al ) and (I4.6bl) are rewritten as 



-^=Pu + f(u), 8(0, X) = J] e ij " kx % n (X), 

n=\ 

OAn 

dT 



= Q n A„ + R n (A), A n (0, X) = v n {X), (n = 1, • • • , N). 



Integrating them yields 



(Y e ijnkx/v Vn) + e PiT - s) f(u(s))ds, 

id J ° 

" T v n + I e®"' 
Jo 



R n (A(s))ds, (n = !,••• ,N), 



(4.18a) 
(4.18b) 

(4.19a) 
(4.19b) 
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which have mild solutions in B r for < T < T . For the function space B' written in the 
X- variable, we introduce the norm || • ||, 7 by 

IMI„= max mp{r, M \d a ip(X)\ }. (4.20) 

0<\a\<r XeRd 

If we put (p{X) := (p(X/rf) = ip(x) for a given <p(x), it is easy to see that \\(p\\ = \\<p\\ n , where 
\\ip\\ represents the standard norm on B r . In the present notation, Prop. 3. 2 is restated as 
follows: there exist to, C\ > such that the inequality 

\\e pT (e ij " kx/r, v n ) - e ijnO>Tls + ij n kxin e Q n T Vn j ^ < ^ j ^ ^ (4 2 \) 

holds for T > st and for each n = 1, • • • ,N. Let us estimate w - ^ =l A n e ijn0Jt+ ' j " kx by 
using the norm || • ||^. 

Step 2. Gronwall inequality. It follows from Eq.(4.18) that 

N N N 

u _ ^ e ij„ojT/£+ij„kX/l — e^ T {e^ nkX ^ ] Vn) — g^*^ '/e+ijnkX/v g &nT ^ 

n=l n=l n=\ 

+ f e PiT ~ s) f{u{s))ds - V C ^I^^^^RjAisMs 
Jo ^ Jo 

= F(T) + ^ e p(J - s \f(u(s)) - f(^n=i A(s)e ijn,JjslE+ij " kx/r >))ds, 

where 

N N 
F(T) = ^ T {e ikkXh V n ) - V e Hn«>Tls+ij n kXlr, e Q n T Vn 

n=l n=\ 

/ ( ^ s) /(Z!-i A(^ !7 «^ /s+l7 "* x/ ^s - V e ,7 "" r/£+,7 "^ /, 'e 4(r -' s) J R„(A(5))rf5. 

° id Jo 

Because of the existence theorem of mild solutions, there exists a positive constant D\, 
which is independent of s, such that 

N 

< D u \\J]MT)e ij ^ T/E+ij " kx % < D u \\e pT \\ J} < D, 

71=1 

holdforO <T <T . LetL > be a Lipschitz constant of / in the ball {<p e B r \ \\<p\\ v < D { ). 
Then, we obtain 

\\u - V A n e ij "" T,E+ij " kx % < \\F(T)\\ V + I DiL\\u(s) - A n (s)e ij "^ /£+ij " kX %ds. 

id J° 

Gronwall inequality yields 



YA t /^ T/E+i J" kx % < \\F(T)\\^ + D X L f e^^WFis^ds. 
7Ti J° 
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To estimate F(T), we rewrite it with the aid of Eq. (l4.2l ) and C jn = R n as 

N N 

F{T) = Yj e pT (e ij " kx/r} v n ) - ^ e ij,^T/ E+ ij„kx/i e Q„T Vn 

n=\ n=\ 
N r T 

+ ^ I ^ j/cr-^Att/i _ e i ^-'V wi i tt /^-'))/{ ) ,(4( J ))J i 

+ J]Hj(T,X), (4.22) 

where i? ; is defined by 

Hj(T,X)= f e p{T - s) e ij " sle+mir, Cj(A(s))ds. (4.23) 
Jo 

Step 3. estimate ofHj. Let C([0, T ]; 5 r ) be a Banach space of functions g(T,X) on 
[0, To] x R rf such that T i-» g(T, ■) e B r is continuous. The norm is defined by 

||g|| c o, := max \\g(T, (4.24) 

Te[0,T ] 

Let C([0, T ]; 5 r ) w be the product space with the norm defined by ||g|| c o.,- = max ||^„|| c o.,- 

l<n<JV 

for g = (gi, • • • , g N ). Due to the existence theorem, a mild solution A(T, X) of (I4.19bl) is 
included in C([0, T Q ];B r ) N . The next lemma will be used several times. 

Lemma 4.4. Suppose / : B r — > 5'' is C 1 . There exists a function /z : [0, T ] x R d x 
C([0, r ]; ST -> C([0, T ]; 5'), which is bounded as 77 -» 0, such that 

£ tf/r, X) = T]h(T, X, A(T, X)). (4.25) 

Further, if / : B> -» 5' is C 2 , /z(r, X, A) is Lipschitz continuous in A e C([0, T ] ; 5 r ) w . 
Proof. In the ^-coordinate, we have 

/(r-V^/ £+ ^C/A>,X)) 

= [^'^^CjiMs^x + riy)) fe'^e PmT ~ s)/£ dyd^. (4.26) 

(2tt) J J 

There exists a number < 6 < 1 such that 

e p( - T - s) e ij0Js/E+ijkX ^Cj(A(s,X)) 

= -J— [e ija)sls+mx+y) C,{A{s,r}x)) f e~^e PmT - s)/e dyd^ 
(2n) d J J 

+ TTtA e ij » s/s+iMx+y) V ^- C 7 (A(5,v)) • y t [ e^e^^dyd^. 
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Since §e~' y ^e p ^ ){T ~ s)le d£; is rapidly decreasing in y, the right hand side above exists. We 

denote the first term and the second term above by I\ and I 2 , respectively; Hj = Chds + 

I 2 ds. At first, we consider I 2 . Since / is a C 1 function on B' and [0, T ] is a finite 
interval, / is regarded as a C 1 function on C([0, T ]; B r ). Since the derivatives dC/s 
are Fourier coefficients of df, the series Yijtj h converges and there exists a function 
h 2 (T,X,A) such that 



r-T 

V I 2 ds = r 1 h 2 (T,X,A). 
"77 Jo 



From the definition, we verify that h 2 is a mapping from [0, T ] x R rf x C([0, T ];B'~) N into 
C([0, T ];B°) (we will show later that this is a mapping into C([0, To];B r )). Furthermore, 
if / is C 2 so that C/s are C 2 , then /? 2 (7/, X,A) is C 1 in A e C([0, r ]; fi r )^ (in particular, 
Lipschitz continuous). 

Next, let us calculate the first term Ii . Note that the equality 



e iiky e~ iy te PmT - s)l£ dyd% = (2n) d e Pm(J ~ s)lE (4.27) 
holds, which can be proved by the same way as Lemma 3.4. Thus we obtain 

fhds = e ijkx e pmTlE f Cj(A(s, V x))e^- pm)s/s ds. 
Jo Jo 

If A is differentiable in s (i.e. when the initial condition is included in the domain of 
<3„X then integration by parts proves that the above quantity is of 0{e). When A is not 
differentiable, we need further analysis. 

Let J' be the set of integers j such that j £ J and Re[P(ijk)] = 0. Due to the assump- 
tion (D3), J' is a finite set. Put ijco - P(ijk) = pj + iqj with pj, qj e R. When j £ J U /' 
(i.e. pj ± 0), the mean value theorem proves that there exists < Tj < T such that 

f Cj(A(s, r]x))e {p J +iq J )s/£ ds = Cj(.A(Tj,T]x))e iq ^ ,£ [ e p ' sl£ ds 
Jo Jo 

= -C ; (A(T ; -,77x))^^ /e (^ r/£ -l). 

Pj 

Since C/s are Fourier coefficients of a C 1 function / and since pj — » oo as \ j\ —* oo, the 
following series 



^ e ijkx e PmTlE —Cj{A{T ] ,T 1 x))e iq ^ l %e^ Tls - \) 

s ^ e iikx e nm[P{imTl£ —C j {A{T j ,T]x))e iq ' T ' l£ {l -e~ p ' T> 
jtjur Pj 
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converges, and there exists a function hi(T,X,A) from [0, T ] x R d x C([0, T ];B r ) N into 
C([0, T ]; 5°), which is C 1 in A, such that 



r 7 

si 



/ids = e^(r,X,A). 



To estimate the case j e J', we need the next lemma. 

Lemma 4.5. For any constants c, t\ > that are independent of s, a mild solution of 
(l4~T9bl) satisfies \\A n (T + cs) - A n {T)\\ n ~ 0{r]) for eh < T < T and n = 1, • • • ,N. 

Proof. Recall x = X/rj, t = T/s. In the (?, x)-coordinates, linear semigroups satisfy 

e &(T + cs) Vn _ e Q„T Vn = C Vn(r]X + j f e -«tf (g Q„(*Xf +c ) _ e QM)t )dyd ^ 

(2n) d J J 

Since 

J v n (^) J e-^e^'dydt = v n (r]x) 
for any t, there exists < 6 < 1 such that 

e Q^cs) Vn _ e Q n T Vn = V Cy + ^ C e -iyt e QM)t {e QMY _ Y)dyd%, 

(2n) d J dxi J 

which is of order 0{rf) uniformly in x and t\ < t. Next, the derivatives satisfy 

^ ( e& ' <T+CE)v » - e& ' (T)v ») = (2^/ v '^ x + W) J [^^-^^^(e 6 "^ - l)dfydf . 
By the same calculation as above, it turns out that 

\\e^ T+CE) v n -e^ T v n \\ n ~0{ri) (4.28) 
for t { < t; that is, for et { < T. Then, Eq. (l4.19bl) yields 

XT+CE 
e Q " (T+C£ - s) R n {A{s))ds 

+ (e Q " (T+CE ~ s) - e^ iT - s) )R n (A(s))ds + ( e Q*V«*-*) _ e ^ (T ~ s) )R n {A{s))ds. 

Jo Jr-e/i 

Using (14.281) . we obtain the lemma. ■ 
Suppose j e J', so that ijto - P(ijk) = iqy, 

f I l ds = e ijkx e P{ijk)T/£ f Cj{A{s,r]x))e iq ' sl£ ds. 
Jo Jo 
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Since the set J consists of all integers satisfying P(ijk) = ijco (the assumptions (D3)), 
qj 4- 0. Put 7 3 j = Cj(A(s, j]x))e' q i s/s ds. Changing the variable s i-> s + en/qj yields 



Cj(A{s-en/qj,Tix))e^ s/g d5. 

en/qj 



Hence, we obtain 

-•En/qj fT+en/qj 



X sn/qj rl +SK/qj 

Cj(A(s, r]x))e iq J s/E ds - J Cj(A(s - sn/qj, r]x))e iq ' sls ds 

+ f (Cj(A(s, rix)) - Cj(A(s - en/qj, rjx))) e iq > s ' e ds 

J en/qj 

r-n/qj ^ nT/s+n/qj 

= e I Cj(A(ss,Tjx))e iqjS ds - e I Cj(A(es - en/qj,r]x))e iqjS ds 

Jo Jt/e 

C T Cj(A(s, rjx)) - Cj(A(s - en/ qj, rjx)) . . 
+ 7]- — e q,l ds. 

Jen/q: ?! 



Lemma 4.5 shows that 

. f Cj(A(s, rjx)) - Cj(A(s - en/qj, r\x)) 

h 3 ,M) := — 

rj 

defines a function from C([0, T ]; B r ) N into C([0, T ]; B r ), which is bounded as r] — > 
0. Hence, there exists a function h 3J (T,X,A) from [0,T ] x R d x C([0,T ];B r ) N into 
C([0,r ];5°) such that 



J hds = rjh 3J (T,X,A). 



Further, if / and Cj are C 2 , h^j{T,X,A) is C 1 in A. Therefore, putting h = h 2 + rj M l h\ + 
lijej' Kj P^es Eq.(g23 satisfying h(T,X A) e C([0, T ]; 5°). 
Let us estimate the derivative. Eq. (l4.261 ) yields 

= fjM'+W^CjiAfaTix + Tiy)) f(^Te^e PmT - s)/£ dyd^ (4.29) 

(2^) rf J J 

Repeating the same argument, it turns out that /z(7\ X,A) 6 C([0,T ]; B r ) if A is in 
C([0, r ]; 5')^. This completes the proof of the lemma. ■ 

Because of this lemma, there exists a positive number D 2 such that 

Wj^HjCT^KriDi (4.30) 
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holds for < T < T . 

Step 4. estimate of F(T). By using Eq. (l4.21l) and (14.301) . we can show there exists a 
positive constant Z) 3 such that ||F(r)||^ < 77D3 for sto < T < T . Therefore, we obtain 



n = \ 

< 77D3 + AL e DlLiT - s) T]D 3 ds + e DiL(T - s) \\F(s)\\ v ds ~ 0(Tf). 

Jetn Jo 



„0)T/E+ij„kX/riu 

n=l 



'sto 

for sto < T < Tq- Changing to the (t, Jc)-coordinate proves Theorem 4.2. ■ 

Proof of Thm.4.3. Let us consider the systems (4.15). Recall that in this situation, 
Cj(Ai, ■ ■ ■ ,Ajv,A_i, • • • , A_jv) is defined by Eq. (l4.11D . S n (A) is defined by 

S„(Ai,--- ,A N ) =R n (A u --- ,A N ,Ai, ■ ■ ■ ,A N ) = C jn (A u --- ,A N ,A U --- ,A N ) 

for n = 1 , • • • , N. Again we assume D = d and use the same notation as the previous 
proof. A mild solution of (14. 1 6al ) written in the (r,X)-coordinate satisfies 

u = e PT {u{0))+ \ e PiT - s) f(u(s))ds, (4.31) 



X 

with the initial condition u(0, X) = u(0). Let us consider the system of 2Af-integral equa- 
tions of w+ := (w\ , ■ ■ ■ , wjv) and w_ := • • • , W- N ) of the form 

e ij»"T/s+ijnkxh Wn _ e ^( e iJnkx/v Wn (Q^ + f e p ^ T ~^R n (w + (s),W-(s))e ijnaJ ^ s+iJnkx/r, ds 

Jo 

+ f e p( - T -%(w + ( S \w-(s))e ij " s/e+ijkx/ «ds, (4.32) 

e -ij H6 >Tj*-ij H Wn w _ n = e pT (e- ij " kxJrj W- n (0)) + f e PiT - s) R^ n {w + {s), w^s))e~ ij " ws/£ - ij " kx/ ids 

Jo 

f / (r - i) C_/w + (5),^(5))e- i > l/£ -'™J5, (4.33) 

which has a unique solution vv„ e C([0,T ]; B r ) satisfying w„(0, X) = w n (0) e B r . This 
yields 



N 



e ij„uT/£+ij„kX/r} w 



n=-/V 

N ^ „j 

/ r ( ^ e'^ x/ V ;i (0)) + ^ I Jv-*) C j( w+ (s),w-(s))J* os/ ' l * 4Jkx, '>ds 

n=-JV y'=-oo ^° 

/ r (Y e ij " kx/T, w n (0)) + / (r - s) /(Z^=-ivW„(5)e 0>i/e+ ^ fcz/ ' 7 M5, 
•f-jl Jo 



n=-N 
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where we used the abbreviation J^=-n := ^n=-\ + 2^=i- Recall also that j_„ = -j n . 
This means that u = Z„=-n e ikoiT,E ^ kX ^w n when u(0) = Zn=~N e ij " kX/T, w n (0). 
By using (D4), we can show the equalities 



and 



Cj(A u --- ,Ajv,A_i,--- ,A_ N ) = C_ 7 (A_i,--- ,A_ N ,A U --- ,A N ). 

Therefore, {w n = w~„ e R}^ =1 is the invariant set of the 2iV-equations. By putting w n 
W- n e R and w := w+ = W-, the system is reduced to Af-equations of the form 

e ij n ojT/£+ij„kX/ri w ^ _ ^PT ^ij^X/rj 



w„(0)) + f e p(T - s) S n (w(s))e ijntos/E+ii " kX,T, ds 
Jo 

i r T 

+— Yj J e ?,(r "' s) Cy(w(^))e , >' s/e+, '^ / "J5, (4.34) 

for n = 1, • • • ,N, where Cj(w) := Cy(w, w). 

Suppose that there exists (p = (fa, • ■ • , (p N ) e R w satisfying the assumptions of Thm.4.3. 
Without loss of generality, we assume that the Jacobi matrix of (5 1, • • • , S N ) at A = <p is 
diagonal. Thus we put 

Q CI 

S n (fa,- ■■, ( f> N ) = 0, jr^ifa, • • • , <p N ) = -J3„ ■ 6 n , m , Re\J3 n ] > 0, (4.35) 

ioxn,m= l,--- ,N. Putw„ = </>„ + r]W n . Due to the assumption of Thm.4.3, S n is C 2 hav- 
ing the Lipschitz continuous second derivatives. Hence, there is a Lipschitz continuous 
function S„( • ,77) : C([0, T ]; B'") N —> C([0, T ]; B r ), which is bounded as 77 — » 0, such that 

+ 77W) = -77^ W n + r] 2 S n (W, 77), W = (Wu ■■■ , W n ). (4.36) 
We denote S n (W, 77) by S„(W) for simplicity. Note that Eq. (l4.27l) gives 

e pT (e' kkx ^(p„) = e P(i ^" k)TlE e^" kxlrt (t) n = e^ nuTl£+i ^" kxlri (p n 
Thus Eq. (l4.34l) is rewritten as 

e ij„o J T/ £+ ij„kxh Wn _ e pT (e ihkxh W n (0))- B n f e p{T - s) W n (s)e ij " tos/£+ij " kx/ ' 1 ds 

Jo 

+ 77 f e p(T ' s) S n (W(s))e'^ s/£+ij " kx/,1 ds + Y f e p(T ~ s) C U + rjW(s))e ijojs/s+ijkx/r >ds. 
Jo 277N j£ Jo 

Remark that the second term in the right hand side is linear in W n . Therefore, we can 
show that this equation is rewritten as 

e ij n coT/ s+ ij n kx/n Wn = e {p ^" )T (e ij " kx/r, W n (0)) + ri f e (p ^" )(T ' s) S n (W(s))e u "^ /£+ij " kx/ri ds 

Jo 

+^A?y f e (p - MT - s) Cj((P + r]W(s))e ijojs/s+ijkxJrj ds. (4.37) 



2riN^ J,, 
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Motivated by this equation, let us consider the system of integral equations of the form 

%J -co 

V f e V>-Pn)(T- s ) e ijcos/ s+ ijkx/ vc / . + ^(s^ds. (4.38) 

We will show later that a solution W* = (W*, • • • , W^) of this system satisfies Eq. (l4.37l) 
with a suitable initial condition W,,(0) = W,,(0, X). 

To prove the existence of a periodic solution, let C p be the set of functions i//(T, X) in 
C([0,T ];B r ) such that 

2ns 

ifr(T + ,X U --- ,X d ) = iff(T,X u --- ,X d ), whenw^O, ,.~ Q . 

ifj(T, Xi, ■ ■ ■ , X d ) is constant in T when co = 0, 

and 

<ff(T, X u ■ ■ ■ , Xj-u Xj + X J+U ■■■ ,X d ) = if/(T, X { , • • • , Xi), when kj * 0, 

if/(T, X\ , ■ ■ ■ , X d ) is constant in Xj when kj = 0, 

(4.40) 

for j = 1, • • • , d. By the norm \\tf/\\c p '■= max reR \\if/(T, • )||^, C' p becomes a Banach space, 
which is a closed subspace of C([0, T ];B r ). Let (C p ) N be the product space with the norm 
\ Cp = max^jv \\ifr„\\ C r p for i/f = (ft u ■ ■ ■ , iff N ). Define mappings Cl Un and Q 2 ,n to be 



(Q 2 ,„W)(r,X) = — J] f e~ 0>r/e ~' 7 "* x/ V^" )( ^ 
for « = 1, • • • ,N. 

Lemma 4.6. and £2 2 ,« are mappings from (C p ) N into C£. 

Proof. By using the expression (13.61) of the semigroup, we can show that if W e (C p ) N , 
then there exists a positive constant Di such that 



^1 



r . e -ijnkX/V e P(T- S ) e ij n kX/ V§nms! X)) 



dX a 

for X 6 R rf , -oo < s < r and |ar| = 1 , • • • , r. Thus we obtain 



wn^wwci < n ■ sup f D ie - MT - S) ds. 

TeR J-co 



Since Re[/3„] > 0, the right hand side above exists. The periodicity conditions (I4.39I) . (I4.4()I) 
immediately follow from the definition. The proof for Q>2,n is done in the same way. ■ 
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Lemma 4.7. Fix a positive number 6 and let D = {W e (C p ) N 1 1| W\\ Cp < S} be a closed 
ball in (Cp^. If 77 = 77(5) > is sufficiently small, • • • , Qi,jv) and ■ • * , ^2jv) 
are contraction mappings on D. 

Proof. It follows from the proof of Lemma 4.6 that ||Oi, n W||c^ is of order 0(rj). Thus it 
is easy to verify that (Qi,i, • • • , Qi,jv) is a contraction mapping on T> if 77 is sufficiently 
small. Next, by the same way as the proof of Lemma 4.4, we can show that there exists a 



function h n : R x R d x (C r n ) N -> C r n such that 



(n 2>n W)(T,X) = -L- • r ] h n (T,X,cf> + T}W(T,X)), (4.41) 

where h n (T,X, ■ ) is Lipschitz continuous. Therefore, there exists L n > such that 

P 2 ,„W - fi^VlIc; < ^\\r]W - tjVWc; ~ 0(tj). 

This proves that (Q 2 ,i, ■ • * , is contraction if 77 is sufficiently small. ■ 

Due to this lemma, the system (14.381) of integral equations has a unique solution W* = 
(W[ , • • • , W* N ) in (Cp w , which is periodic in X and T. 

Lemma 4.8. The solution W* is a solution of (14.371) satisfying the initial condition 

f° - 

jJ*xx/vW n (0) = ^ e^- M 'J i ' as/ *+ ii »* x, ' l S n Wr(sy)ds 

*J —00 

1 r° 

V e HP-fi«)s e ija>s/e+ijkXM C y + rjW*(s))ds. (4.42) 



Proof. This follows from the substitution of (HT4"2l into (14371 ■ 

Now we have proved that the system (14.341) has a solution w n (T,X) = <f> n + r]W*(T,X) 
satisfying W* e C p . Therefore, the equation (14.311 ) has a solution 



u(T, X) = y {(f> n + r]W*(T, X)) ■ (eijn<oTis+ij n kXM + e -ij n ajT/ s -ij n kx/ V y (4 43) 

n=l 

Changing to the (t, x) -coordinate yields a mild solution of (I4.16al) of the form 

N 

u(t, x) = J^ (<p n + r}W* n (st, Tjx)) ■ { e iinU)t+ijnkx + e ~ ij " a "- ij " kx l (4.44) 



which proves the first part of Thm.4.3. 

Finally, let us prove the stability part of Thm.4.3. Let W*(T, X) be the periodic solution 
of (14371) . There exists a positive constant D x > 1 such that \\ e -^hMh e PT e ihMh^ < ^ 
for any T > and n = 1, • • • , N. Fix a positive number M and put 5 = M/(2Di). 
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Let D = {W e (B r ) N \ \\W - W*(0, • )||^ < 6} be a neighborhood of the periodic solution 
W*(0,X) at T = 0. Due to the existence theorem of mild solutions, there exists To > 
such that when W(0) e D, then Eq. (l437T) has a solution W(T) in (5^ for < T < T . 
We define a time To map as 

T : £> -> (B r f, W(0) h> W(r ). (4.45) 

Lemma 4.9. If r\ > is sufficiently small, To can be taken so that T is a mapping on D. 

Proof. Define Q 3 „ and £\„ to be 

Jo 

(Q 4n W)(T,X) = — Y f e-^^^-^^^e^^^-^e^^^^^C/^+nWC^.Z))^, 

for n = 1, • • • , N. Eq. (l4.37l) gives 

W„ - = e- iJ " uT '^ i ^ , *eP^ T j J ^ / ''(Wn(0) - W*(0)) 

+r 1 (^ n w - a Xn w*) + -(O^w - n 4>n w*). 

By the same way as the proof of Lemma 4.4, we can prove that there exist Lipschitz 
continuous functions h n : [0, T ] xR''xC([0, T ], B r ) N C([0,T ],B r ) such that Q 4 ,,rW = 
r]h n (T,X, (p + T]W). This provides 

W n - W* = _ w*(0)) 

+T](Q Xn W - Q Xn W*) + hn(T, X,(f> + TjW) - hn(T, X, (f> + T}W*). 

Since ||W(0) - ^(0)11, < 6 = M/(2D{) for W(0) e £>, we can assume that T is chosen 
so that \\W(T) - W^T)^ < M for < T < T . Put B := min l£n < N fi n . Since §„ and h n are 
locally Lipschitz continuous, there exist positive constants D 2 , D 3 such that 

(? T \\W{T) - W*(T)\\,, < M/2 + r]D 2 [ (P s \\W(s) - W*(s)\\ v ds + rjD^^WW - W*\\co, 

Jo 

for < T < To (see Eq. (l4T24l) for the definition of the norm || • || c o, on C([0, T Q ];B r ) N ). 
The Gronwall inequality gives 

e^\\W(T) - W*(T)\\ n 
< M/2 + TjD 3 e^ T \\W - W*\\ c o,,- + r]D 2 f e^ 2 ^ (m/2 + rjD^W - W*\\ c o, r ) ds 

= Me^ T /2 + -H£^-e^\\W - W*\\ c o, r - !J^l e ^ T \\W - W*\\ c o, r 
B- n D 2 MC B-t]D 2 

for < T < Tq. Hence, we obtain 

\\W(T) - W*(T)\L < Me illD2 - p)T 12 + \\W - W*\\ c o, r . 

S-T]D 2 
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Using the standard existence theorem, we can verify that a solution W n (T, X) of (14.371) is 
bounded as rj — > 0. Thus, \\W - W*\\ c o.r is bounded as rj — » 0. Therefore, if rj is sufficiently 
small, we obtain ||W(r) - W*(r)||^ < M/2. This implies that T can be taken arbitrarily 
large and \\W(T) - W*(T)% < M holds for any T > 0. Hence, if T is sufficiently large, 
we obtain ||W(7o) - W*(T )\\ ri ~ 0(rj). Since W* is 2^e/a»-periodic in T, we can choose 
T so that \\W(T Q ) - W*(0)\\ v ~ 0{rj), which proves W(T Q ) 6 £>. ■ 

Lemma 4.10. If rj > is sufficiently small and T is sufficiently large, T is a contraction 
mapping on £). 

Proof. Let W(r) and V(T) be two solutions of (14371) with the initial conditions W(0), V(0) e 
£), respectively. By the same calculation as above, we can show the inequality 

\\W(T) - V(T)\\ n < Di\\W(0) - V(0)\\ n e (flIh - p)T + rjD 4 \\W - V\\<*, 

where we put D 4 = D 3 j3/(j3 - J]D 2 ). Since \\W - V\\ c o, = max <r<r W(T) - V(T)%, we 
obtain 

\\W - V\\co, < Di||W(0) - V(0)||, + ?]D 4 \\W - V\\ c o, r . 
Substituting this into the above inequality provides 

\\W(T) - V(T)% < \D x e^-® T + Y^^j W(0) - 7(0)11,, 

which proves the lemma. ■ 

Take T = 2ns /to ■ I so that Lemmas 4.9 and 4.10 hold, where / is a sufficiently 
large integer. Due to the lemma, there exists a unique function W**(X) e D such that 
any solutions W(T,X) in D satisfy W(nT , •) -> W**(-) as n -> oo. If W** * W\ then 
W*(nT , ■ ) converges to W** as n — > oo. This contradicts with the fact that W*(nT , ■ ) = 
W*(0, ■ ) is independent of n. Hence, W** = W* and W(T,X) in D converges to W* as 
T —* oo. Since a mild solution of (I4.16al) is written as 

N 

u(t, x) = J] (<f>n + rjW n (st, 77*)) ■ + g"^^**), (4.46) 

the proof of Thm.4.3 is completed. ■ 
4.2 Higher dimensional case 

Suppose « = («!,-••, w m ) G C m and x = (xi, • • • , xj) e R f/ . For fixed 1 < D < d, we 
use the same notation x = (xi,x 2 ), a = (/3, y) as in Sec.4.1. Let {P,y(x)}™. =1 be the set 
of polynomials of x. The mx m matrix P(x), the differential operator V and eigenvalues 
• • • , /l m (^) are defined in the same way as Sec. 3. 2. We suppose for simplicity that 
only A\(£) contributes to the center subspace of V (see (El) below). Extending to more 
general situations is not difficult (see Remark 3.7 and Example 4.12). 
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(EO) The matrix P(ig) is diagonalizable for any £ e R d . 

(El) Re[Ai(&] < and Re[i/(£)] < for any £ 6 R d and j = 2, ■ ■ ■ , m. 

(E2) There exist co e R, k e R f/ ((a>, k) 4- (0, 0)), a finite set of integers J = iji, ■ ■ • , j N } 

and {Mi, • • • , M N ) such that 

Mjnk) = ijnio, (n = 1, •• • ,N), 



as* 



(j n k) = 0, for any B such that \B\ = 1, • • • , M n - 1, (n = 1, • • • ,iV), 
0'„/c) ^ 0, for some B n such that \B n \ = M„, (n = 1, • • • ,N). 



The set J consists of all integers satisfying A\(jk) = ijco. 
(E3) For n = 1, • • • ,N, define Q n {x) and Q n by 

QnW = Qni&uG) = y „ 1 „ A U j n k)(x X lif, Qn = Qn(dy , ■ • ■ , d D , 0, ■ • ■ , 0). 

(4.47) 

Then, both of 'P and <3„ are elliptic in the sense that there exist C\,c 2 > such that 
ReW;(#] < -c 2 |f| 2 (j = l,-- ,m) and Re[Q„(/|i)] < -cz&fCn = l,~-,N) hold for 
Ifl.llil >c l5 where & = (ft,.-- ,£,). 

When m = 1, /li(£) = P0'£), so that the above assumptions and <2„ are reduced to 
those given in Sec.4.1. Let 5' = BC'(R' ; C) and (fi r ) m = B r x • • • x B r a product space. 
The norm on (5 r ) m is defined by \\u\\ = maxi<„< m ||w„||. Note that P is an operator densely 
defined on (B r ) m , while Q„ is an operator densely defined on B r . 

Let w n = (w„ t i, ■ ■ ■ ,w n j n ) be an eigenvector of P(ij n k) associated with the eigen- 
value Ai(j n k) = ij n co for n = 1, • • • ,N. The projection to the eigenspace span{iv„} is 
denoted by n„, and the projection to the eigenspace associated with the other eigenvalues 
AiUnk), • • • , A m (j n k) is denoted by IT ; f = id - H n . Functions e ljlkx Wi, ■■■ , e' jNkx w N span 
the center subspace of P. 

When k = (ki, • • ■ ,kd) 4 0, we can assume without loss of generality that k\ ± and 
ji = 1. For a given function / : (B r ) m -> (B'J n , define a function C 7 : (B r ) N -> (fl'") m by 



Cj(A) = Cj(A u ...,A N ) = ^-J^ f(2j A„e ijnhxi w n )e~ ijklXl dxi. 



(4.48) 



n=\ 

When /c = and <x> ^ 0, we use 

■*2n/co N 



r>zn/co " 

Cj(A) = — f(yA n e ij » M w n )e- ij0it dt (4.49) 
instead of (14.481) . Then, we obtain the expansion 

TV oo 

/(^ A n e ijntM+ijnkx w n ) = Yj C M) eii0>t+iikX - ( 4 - 5 °) 
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Further, define the function R n : (B r ) N -> (B r ) m to be 

R n (A) = U n Cj„ (A), j„ e J, (n = 1, ■ ■ ■ , N), 



(4.51) 



and define R n : (B r ) N — > B' so that R n (A) = R n (A)w n . Let s > be a small parameter. Let 
us consider the two initial value problems: 



du 



N 



— =Pu + sf(u), w(0, x) = 2_j e' ] " kx v n {rixx)w n , 



dt 
dA 



(4.52a) 



n=\ 



- = Q n A n + eR„(A), A„(0, x) = v ;i (^i ), (n =!,-■■ ,N), (4.52b) 



I dt 

where 7/ = e 1/M and M := min{Mi, • • • , M N ). 

Example 4.11. Let us consider the system (I1.5I ). whose perturbation term is given by 

f(u, v) = 



u — u 




(4.53) 



The reduction of the linear part was calculated in Example 3.5, in which it was shown that 
(EO) to (E3) are satisfied with m = d = 2, D = 1, k = (0,c), co = 0, J = {j\ = 1 , j 2 = — 1 } 
and Mi = M 2 = 4. We use the same notation as Example 3.5. The matrix P(z£) at 
= (0, ±c) is given by 



P(+ic) = 



(k + d)/(2d) -1 

1 (k + d)/2 ) " 



Eigenvalues and eigenvectors of this matrix are 



w 



V = 



1 

(k + d)/2 
1 

(k + d)/(2d) 



for Ai(0, +c) = 0, 

f + ^ q-d)(k + d) 

for A 2 (0, ±c) = — < 0. 

2d 



Since both eigenvectors of ii(0, c) and A\(0, -c) are given by w above, we calculate the 
Fourier expansion of f(A\e tcy w + A^e'^w). Then, it turns out that 



d(A) 



Aj - 3A^A 2 




C_i(A) 



A 2 - 3A!A2 




Then, it is easy to show that projections of them are 

|2 



R,(A) = n!d(A) 



Aj-3AfA 2 A 2 -3AiA^ 



1 - d 



-w, R 2 (A) = n 2 C_i(A) 



1 - d 



■w. 



Therefore, the amplitude equation (I4.52bl ) is given by 



' dAx ^. s 



dt 



= QA! + 



{ dt 



1 - d 

s 

1 -d 



(A X -2>A\A 2 ), 



dA l=QA 2 + ^— 1 (A 2 -3A l A\), 



(4.54) 
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where <3 is defined by (13.171) . If we suppose A 2 = A\, the equation (11.231) is obtained. 

Example 4.12. Noting Remark 3.7, to extend (El) to the case that several eigenvalues lie 
on the imaginary axis is straightforward and the amplitude equation is defined in a similar 
manner as above. Let us consider the system (11.71 ). whose perturbation term is the same 
as the above Example. By the same calculation as Example 3.8 and 4.11, we obtain the 
amplitude equation 

—l=QA l + -(A l -3A 2 l A 2 ), 

PI 2 (4.55) 



dt 



QA 2 + -(A 2 - 3A 1 Ai), 



where Q = Dd 2 is obtained in Example 3.8. If we suppose A 2 = A\, the equation (11.81) is 
obtained. 

Theorem 4.13. Suppose (EO) to (E3), / : (B r ) m -» (B r ) m (r > 1) is C 1 and s > is 
sufficiently small. For any {v„}^ =1 c B' , there exist positive numbers C, To and to such that 
mild solutions of the two initial value problems (4.50) satisfy 

N 

\\u(t, x)-Y A n (t, x)e ij>M+ij " kx w n \\ <Ct] = Ce l/M , (4.56) 

n=\ 

for to < t < T /s. 

Further, suppose that 
(E4) For |0| = M n and n = 1, • • • , N, 

Ph.... 



P(i& = P(-ia f(u) = f(u) and -^fUnk) = -^i-Jnk). 

In this case, the set J consists of J = (ji, • • • ,j N } U • • • , -j N \ as is Sec. 4.1. We 
consider the two initial value problems: 

( n N 

^=Pu + sf(u), «(0, x) = J] (e ij " kx + e- ij " kx ) v^rjx^w,,, (4.57a) 
8A„ 



dt 

Q n An + sS n (A), A„(0, x) = v n (Tj Xl ), (n = 1, • • • , AO, (4.57b) 



v dt 

where S n (A) is defined by (14.151) . In the next proposition, B r = BC'(R d ; R) denotes the 
set of real-valued functions. 

Theorem 4.14. Suppose (EO) to (E4) and / : {B r ) m -> (B r ) m (r > 1) is C 2 such that the 
second derivatives are locally Lipschitz continuous. Suppose that there exists a constant 
vector <p = (0!, • • • , <p N ) e R N such that 

(i) S„(0) = 0forn = l,--- ,N, 

(ii) the Jacobi matrix of (S\, ■ • ■ ,S N ) at (f> is diagonalizable and all eigenvalues of the 
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matrix have negative real parts. 

If e > is sufficiently small, Eq. (l4.57al) has a solution of the form 



U p (t, X) = YfanWn + Ttfjt, X, T})) ■ (g^^-** + g^**-***). ( 4 .58) 



The vector- valued functions ift n and w p are bounded as 77 — » and satisfy 

2;r/a>-periodic in t (when a> £ 0), J 2^/fcy-periodic in x 7 - (when fcy ^ 0), 
constant in ? (when co = 0), | constant in xy (when fcy = 0), 

for j = 1, • • • , J. This w p is stable in the following sense: For any n = 1, • • • ,N, there is a 
neighborhood U„ c 5 r of <p n in 5'' such that if v„ £ U n , then a mild solution u of the initial 
value problem (14.57al ) satisfies \\u(t, •) - u p (t, -)|| — » as t — » 00. 

Proof. We suppose D = d for simplicity and use the same notation as the proof of Thm.4.2 
(see Step 1). Two mild solutions of Eq.(4.50) satisfy 



U - V &JhuTls±iinMh yVn = V e PT( e ij n kX/r, VnWn) _ y e ij n a>T/ S+ ij n kX/r,( e Q n T Vn)Wn 
n=\ 11=1 n=i 

e PiT - s) f(u(s))ds - V e ,7 "" r/£+,7 "^ / "e 4(r -' s) J R„(A(5))H'^5 

= F(r) + J e p(T - s) (f(u(s)) - f(Zn=i A{s)e iina>slE+ii " kxlll w H ))ds, 
<J 



where 

N N 



F(T) = ^ e pT (e ij " kx/T! v n w n ) - ^ e ik(aTl ^ ikkx, \e &nT v n )w n 

n=\ n=\ 

e ^- S )f^N =i ^ s yUo i sis^ij n kXM Wn)ds _ V e ij ^ T/£+ij " kx,,1 e aAT - s) R n (A(s))w n d S . 
Trf Jo 



As before, Gronwall inequality yields 



£=1 Jo 



||F(5)||„J5 
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for some constants D u L > 0. By using definitions of Cj and R n , we rewrite F(T) as 

N N 

F(T) = e pT {e ijnkXlT, v n w n ) - ^ e ii * a,T/ ** iJ "* x, ' l (e** T vJw n 



n=l n=l 
N n T 



+ J]hj(t,x) 

+ y e ij ^ sls e PiT - s) e ikkXl ^C jn (ks))ds, (4.59) 
where Hj is defined by 

Hj(T,X)= f e? (T - s) e ij0,s/e+ijkx/T! Cj(A(s))ds. (4.60) 
Jo 

Now we have arrived at the same situation as (14.221 ) except for the last term. 

Lemma 4.15. There exists a function g n : [0, 7b] x R d x C([0, T ];B r ) N -> C([0, T ];B r ) m 
such that 

f &™/*/W-W*x/>^ CjH( A (s))d5 = ng n (T,X,A). (4.61) 
Jo 

g(7\ X, A) is Lipschitz continuous in A e C([0, T ]; 5 r ) w and bounded as rj — » 0. 

If this lemma is true, the rest of the proofs of Thm.4.13, 4.14 are completely the same 
as those of Thm.4.2, 4.3. 

Proof. We calculate e p(T - s) e^ kx/ mf;Cj n (A(s)) as 

e p(T-s) e u t M/ m ± c {A{s)) = jl ff^v^^^^^^^ 

(2tt) J J 

Jj„kX/rj n n 

= W J J e ^^ +iMXr_i)/ X C,,(A(,,X +3 0)^. 

Let us suppose that the coordinate of u = (u\, • ■ • , u m ) is defined so that the m x m matrix 
P(ij„k) is diagonal. Define a matrix S(£) such that Eq. (l3.22l) holds. By the assumption, 
S (j n k) is the identity matrix. Then, we obtain 

e p(T - s) e ij '' kxl ^ n C h {A{s)) 

JjnkX/rj r* n 

By expanding S(j n k + r]^) ±x , it turns out that there is a function G n {T,X,A), which is 
Lipschitz continuous in A e C([0, T ]; B r ) N and bounded as r\ — » 0, such that 

Jj„kx/v rr 

= W J J ^ C ^( A ^' x + + ^„(t - x, a). 
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Let be the l-th component of the vector Cj H . Due to the definition of JT^, the first com- 
ponent of e Mjnk+v&(T-s)/ SI[ L C ^ is zero? and the i_ th coni p 0nent i s given by e *iUnk+vf)(T-s)/e C (j) 
for / = 2, • • • , m. To prove the lemma, it is sufficient to estimate 



T 

ij„ws/s 



Because of (El) and (E3), there exists a positive number fi such that 

Rc[A,(j n k + nO + P\ < 0, Re[Ai(j n k + tf) + fi] ~ 0(-\{\ 2 ) (4.62) 
as |£| — > oo for any I = 2, ■ ■ • , m. Then, we obtain 

. = f e ij n cos/s e -{i(T- s )/ SKni ( T _ s X ,A)ds, 
Jo 

where ^T„/ is Lipschitz continuous in A e C([0,T ]; B r ) N and bounded as 77 — » 0. The 
mean value theorem proves that there exists < r < T such that 

I n>l = e ij "^ s K„ tl (T - t,X,A(t,X)) C e-^ T ' s)l£ ds, 

Jo 

which is of order 0(e). Hence, putting (0, /„,2, • • ■ , 4,iv) = £^«(7\ A) and 



^g„(r, X, A) = s——I n (T, X,A) + t] [ e ij "^ /E G n (T - s, X, A)ds 
K^n) Jo 

proves the lemma. 



Now the function F(T) in Eq. (l4.59l) is estimated with the aid of Prop. 3. 6, Lemma 
4.4 and Lemma 4.15 to show ||F(r)||^ ~ 0(rj). Then, the Gronwall inequality proves 
Thm.4.13. A proof of Thm.4.14 is also done in the same way as that of Thm.4.3. ■ 
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